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Abstract 

We consider the shear flow of well-ahgned one-component smectic phases, 
such as thermotropic smectics and lamellar diblock copolymers, below the 
critical region. We show that, as a result of thermal fluctuations of the layers, 
parallel (c) alignment is generically unstable and perpendicular (a) alignment 
is stable against long-wavelength undulations. We also find, surprisingly, that 
both a and c are stable for a narrow window of values for the anisotropic 
viscosity. 

PACS numbers:61.30.Gd, 83.70Jr, 83.70.Hq, 47.20.Hw 



1 



Typeset using REVT^ 



In the presence of simple shear flow, smectic phases exhibit a surprising degree of com- 
plexity. As shear rate and temperature are changed, a variety of transitions in orientation 
and morphology have been observed. Although one might expect the liquid layers to simply 
slide over each other with their normals parallel to the shear gradient, the c orientation (Fig. 
|Tp, they often orient with their layer normals pointing in the vorticity (or neutral) direction, 
the a orientation (Fig. ||). This behavior is seen under some conditions in both thermotropic 
smectics and lamellar phases of diblock copolymers f^-^. For thermotropics near the 
nematic-smectic transition, it has been shown that, as a result of nematic fluctuations, a 
alignment is favored over c 0. In diblock copolymers, it has been shown that the wave vector 
dependence of the quartic coupling in the Hamiltonian describing the order- disorder transi- 



tion favors a [T^ . In this paper we consider the steady shear flow of generic one-component 
smectic phases at temperatures well below the critical regime. We show that well-aligned 
(i.e. defect free) systems favor the a orientation: the c orientation, as well as orientations 
intermediate between a and c, suffer an instability from long-wavelength undulations. We 
also find, surprisingly, that within a small window of values for the anisotropic viscosity, 
both a and c are stable. 

It is clear that smectic phases will align so that the average flow velocity has no com- 
ponent along the layer normals. Otherwise, the layers will be forced to deviate from their 
preferred spacing, which is energetically costly. For perfectly flat layers, both a and c orienta- 
tions permit steady shear flows with the layer displacement unperturbed from its equilibrium 
value. However, thermal fluctuations of the layers are convected differently in the two cases. 
As we will see below, convection leads to a greater suppression of thermal fluctuations in 
the c orientation than in the a orientation. Hence, if we adopt a naive picture in which 
the steady-state dynamics is determined by maximizing layer fluctuations (free energy min- 
imization), the a orientation will be obtained in steady state. Since entropy maximization 
arguments in non-equilibrium systems are often suspect, we also compute the dynamic re- 
sponse function for small perturbations away from an aligned state and show that the c 
orientation is indeed unstable towards rotation to a. 



Well below the ordering transition, where the amplitude of the order parameter (concen- 
tration variations for diblocks or density variations for thermotropics) is fixed, a well-aligned 
smectic is parametrized by a layer displacement u{r). We take the average layer normals 
parallel to z. The Hamiltonian for layer fluctuations is (in fourier space) |12[ 

'H = \j d\u{-q)u{q)uj{q) (1) 
uj{q) = Bql + Kqi 

where B and K are, respectively, the layer compression and bending moduli and qx_ denotes 
the component of the wave vector q in the (x, ?/)-plane. The penetration length A = \J Kj B 
is typically of order the layer spacing. As a result, modes with = 0, which correspond to 
layer undulations, are much softer than layer compressions, which have g± = 0. 

Consider now the effect of a steady shear with average flow velocity parallel to x 

v{f) = [f ■ h) jx (2) 
n = cos 9 z + sinO y. 

The velocity gradient direction is n; in the a orientation h = y and in the c orientation n = z 
(Fig. |l]). A mode with wave vector q at time t = will be convected by the flow Eq. (H) 
according to 

q{t) = q- h^q^t. (3) 

In the a orientation, the z component of q is unaffected by the shear, qz{t) = qz', & 
mode that is a pure undulation (g^ = 0) at t = remains a pure undulation at later times. 
In the c orientation, however, qz(t) = 'yqxt; undulations with 7^ at t = pick up a 
compressional character (g^ 7^ 0) at later times. We will see below that the lifetime for a 
fluctuation with wave vector g* is given by l/{uj{q)j3{q)) where /5(g) is a kinetic coefficient. 
As noted above, layer compressions are much stiffer than undulations, i.e. for fixed |g|, uj{q) 
is larger for compression modes than for undulations. Therefore, the lifetime of modes with 
g^ 7^ is much shorter in the c orientation than in the a orientation. (Here we assume an 
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isotropic kinetic coefficient for simplicity.) With a greater suppression of modes in case c, 
we expect a corresponding increase in the "free energy" and hence an instabihty towards 
rotation to a. (It is easy to see that orientations intermediate between a and c suffer from a 
similar suppression of modes, though to a lesser extent, and thus the "free energy" should 
monotonically decrease in passing from c to a.) 

In the previous discussion we have invoked free energy minimization in a non-equilibrium 
system. One can often determine the steady-state behavior in such cases by minimization 
of a quasi free energy. However, this usually requires that the equations of motion can be 
written in the relaxational form cf) = 5T /5(f), where is the dynamical variable of interest and 
r(0) is some functional. In our case, one can consider the equation of motion for the angle 9 
parametrizing orientations between a and c (Fig. m) after averaging over layer fluctuations. 
Unfortunately, we have not been able to express the equations of motion in a form that 
would justify a minimization principle. However, the above picture suggests an alternative 
course. A long-wavelength undulation in the y direction locally looks like a change in 9, 
i.e. a tilt of the layers in the neutral direction. Therefore, a local driving force for rotating 
from c to a should appear in the c orientation as an instability towards long-wavelength 
undulations with g || y. 

Below, we ffist compute a quasi free energy to demonstrate the mode suppression de- 
scribed above. We then compute the dynamic response function (to one loop in perturbation 
theory) and show that indeed the c orientation, as well as orientations intermediate between 
a and c, are unstable and the a orientation is stable. 

To describe the hydrodynamics of smectics, we follow the treatment of |1TT|JT^ . We 



neglect inertial terms and assume incompressibility. In this limit, the dynamics of u reduces 
to the relaxational form: 

^..fe)^^0, (4) 



T] and 7]' are viscosities |T3|, and bp is the permeation constant. Although our results are 
independent of bp, we include permeation in order to ensure that subsequent expressions 
converge for large qz- 

In the presence of steady shear, the time derivative in Eq. is replaced by the convective 
derivative d/dt — 'jq^d/dgn, qn = q ■ n. In order to take into account thermal fluctuations, 
we add to the right hand side of Eq. a random noise, ({q, t), with correlations that insure 
for zero shear rate the system relaxes to equilibrium [Q. We thus have 



- iq.-Q^)u{q) + (3{q)^{<l)u{q) = C(g, t) (6) 

(C(g,t)C(-g,0)) = 2fc5T/3(g)5(t), (7) 

where fc^ is Boltzmann's constant and T is temperature. From Eq. (^, we find for the 
equal-time correlation function of u: 

X{q) = {u{q,0)u{-q,0)) = keT f dt2(3{q{t)) 

xexp(-2^°dt'/3(g(t')V(g(0))- (8) 

Similar expressions are found in [^^-[TB|. 

Since the equation of motion Eq. (|) is linear, x determines the probability distribution 
of u. We thus define a "free energy" density by 

^^-^ /d=^glogx(g). (9) 



2 

Eq. d^) is a complicated function of 7 and 6, however we can extract the leading asymptotic 
behavior as 7 — 0, where 7 = 'yf]/B. For thermotropics, with typical values of -B = 
10* erg/cm^ and rj = 1 poise small 7 implies shear rates 7 < 10* s~^, which easily 



encompass the range studied experimentally. For block copolymers, however, B can be much 
smaller and rj much larger by many orders of magnitude. For example, for B = 10^ erg/cm^ 
and rj = 10^ poise we must take 7 < 1 s^^. We find for small 7 

^ = ^1^=0 + (7 cos 9)'/' + ..., (10) 
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where ci ~ 0.25 is a constant determined numerically. T will therefore be minimized for 
Q = tt/2, which is the a orientation. 

While the above analysis is intuitively appealing in its treatment of thermal fluctuations, 
as we have discussed in the introduction, it is not well justified. We therefore turn to a 
computation of the dynamic response function. 

We can model a time dependent disturbance to the system by shifting the noise in Eq. 
(i) Ci^jt) ~* Cil^t) + fil^t)- The response of the system to f{q,t) is given by 

{u{q,t))=f dt'C{q,t,t')f{q{t'-t),t') (11) 
C{q,t,t') = exp J" ^dt" (3{q{t")Mq{t"))^ . (12) 

For all orientations 6, disturbances created by / relax exponentially with a decay rate that 
depends on the convected value of q. In order to search for an instability, we must go beyond 
the linear Eq. (|^). The smectic Hamiltonian with the leading anharmonic corrections is given 

hym 



+ C{d,u + ^{V^u)% (13) 



The last term in Eq. (|I3]) is a counterterm that is used to enforce the condition 

(£)l/^o = 0, (14) 

which ensures that u describes deviations from the average layer spacing. 

Although we have had to include nonlinear corrections in Ti, this does not mean that 
the physics underlying the response function differs from the physics presented in the intro- 
duction and contained in JF, which was only computed for the quadratic Hamiltonian Eq. 
(0). The non-linear terms in Eq. (p!3| ) are not arbitrary: they are determined by requiring 
that rotating the layers (in the absence of shear) does not cost any energy |2y] . Since JF has 
been calculated for arbitrary tilt 6 within the harmonic theory, it in principle contains the 
same information as the anharmonic terms to one loop. 
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The equations of motion are now nonlinear 

(| - + /^(^) ^ = C{q, t) + f{q, t), (15) 

but we are free to choose f{q,t) arbitrarily small so that linear response is still correct. For 
zero shear 7 = 0, the response function, to one loop order, takes the form 

C{q,t,t') = 

exp (-^°_^dt"/5(g(t"))(a;(g(t")) +5a;(g(t")))) • (16) 

The shift Suj{q) corresponds to a renormalization of B and K; rotation invariance forbids 
a term proportional to from appearing [20|. In the presence of shear, however, we can 



no longer invoke rotation invariance and we expect terms quadratic in q^ and qy to appear. 
For sufficiently long times, the qx dependent terms will be small compared with the Bq"^ 
and terms in u{q(t)). (g^ and qy pick up q^ dependence through convection, Eq. (^). 
A contribution to 6uj that is proportional to and negative, on the other hand, will be 
destabilizing. As discussed in the introduction, this is in accord with the intuition that a 
local tilt of the smectic layers in the neutral direction, i.e. a local change in 9, corresponds 
to a mode with q\\y. 

We have computed the response function C{q,t,t') for Eq. (|l^) perturbatively to one- 
loop order. If we assume f{q,t) is nonzero only for q = qytj, then, in the limit of long times 
(low frequencies) , linear response still takes the form of Eq. ( pA]) with the response function 
as in Eq. (|l^). After taking care to maintain Eq. (|T^) , we find, again in the limit of small 7, 

6uj{q) ^ -c,ql (2 + ^)^(^cos^)^/3 + . . . , (17) 

where C2 ~ 1.4 x 10"'^ is determined numerically. For q along y, uj{q) goes as and so for 
sufficiently small qy, 600 » uo. Hence, for orientations with 9 not too close to 7r/2 (i.e. away 
from the a orientation), disturbances with g || y will grow and the smectic will be unstable. 
Here we have implicitly assumed 2 + rj' /rj > 0; we will return to this point below. 
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When 6 is sufficiently close to 7r/2 such that cos 6' ~ 7sin^, subleading terms in Eq. (|Tj 
will be important and can change the sign of Su. To check the stability of the smectic in 
this limit, we have computed the response function for 6 = 7r/2. In this case we find 

5coiq)^Csql^^ + ..., (18) 

where C3 ~ 4.6 x 10^'^. Thus, for the a orientation, the correction 6uj is positive and the 
smectic is stable. 

As was mentioned previously, in concluding that the c orientation is unstable we have 
assumed 2+//'/ 77 > 0. In fact, there is a small window in which the sign may be reversed. The 
five viscosities that enter the constitutive relation between stress and strain rate in a uniaxial 
fiuid are constrained by the requirement of non- negative energy dissipation |T^ . In our case 



T3| , this requires rj > and rj' > —At]. Thus, as 77' is lowered into the range —At] < f]' < —21], 
we find a transition to a regime in which the smectic is stable for all orientations 6. This 
range of viscosities corresponds to small dissipation for uniaxial extensional fiow along z 
compared with the dissipation for shear fiow. Unfortunately, we have no simple physical 
argument for this result. 

We have considered the shear fiow of well-aligned one-component smectic phases outside 
of the critical regime. Very close to the nematic-smectic transition, where nematic fiuctua- 
tions are large, our analysis is inappropriate. However, the appearance of the a orientation 
in this regime has been accounted for in Our work is similarly complementary to that 
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Toll , where the role of amplitude fiuctuations is considered. Also note that our analysis 
requires modification for two-component systems (i.e. lyotropic smectics) where there is an 
additional hydrodynamic variable. 

We have argued that, as a result of convection and the higher energetic cost of compress- 
ing smectic layers compared with bending, there is a greater suppression of fiuctuations in 
the c orientation. With the naive view that the steady-state behavior is determined by 
minimization of a "non-equilibrium free energy", the a orientation will be favored over c. 
To demonstrate the scenario suggested by these arguments, we have computed the dynamic 

8 



response of the system to a long-wavelength perturbation corresponding to a local tilt of the 
layers. We find that over most of the range of allowed values for the anisotropic viscosity, 
the c orientation and orientations intermediate between c and a are indeed unstable and the 
a orientation is stable. Surprisingly, we have also found that there is a window of values 
for the viscosity in which all orientations are stable. Our treatment completely neglects the 
role of defects as well as the possibility of a nonlinear relation between stress and strain 
rate (non- Newtonian behavior), both of which are likely to play an important role in some 
ranges of temperature and shear. While our analysis by no means accounts for the entire 
phase behavior of smectics under shear, we suggest that the mechanism described in this 
paper may account for the prevalence of the a orientation observed in the shear flow of 
one-component smectics. 
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FIGURES 

FIG. 1. Schematic of the c, a, and intermediate orientations. We always take the flow velocity 

to be in the x direction and the layer normals to point in the z direction. The shear plane is in 
the plane of the page. 
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